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Abstract — Supporting massive device transmission is challeng- 
ing in Machine-to-Machine (M2M) communications. Particularly, 
in event-driven M2M communications, a large number of devices 
activate within a short period of time, which in turn causes 
high radio congestions and severe access delay. To address this 
issue, we propose a Fast Adaptive S-ALOHA (FASA) scheme 
for random access control of M2M communication systems with 
bursty traffic. Instead of the observation in a single slot, the 
statistics of consecutive idle and collision slots are used in FASA to 
accelerate the tracking process of network status which is critical 
for optimizing S-ALOHA systems. Using drift analysis, we design 
the FASA scheme such that the estimate of the backlogged devices 
converges fast to the true value. Furthermore, by examining 
the T-slot drifts, we prove that the proposed FASA scheme is 
stable as long as the average arrival rate is smaller than e^, 
in the sense that the Markov chain derived from the scheme 
is geometrically ergodic. Simulation results demonstrate that the 
proposed FASA scheme outperforms traditional additive schemes 
such as PB-ALOHA and achieves near-optimal performance 
in reducing access delay. Moreover, compared to multiplicative 
schemes, FASA shows its robustness under heavy traffic load in 
addition to better delay performance. 

Index Terms — M2M communications, random access control, 
adaptive S-ALOHA, drift analysis, stability analysis. 



I. Introduction 

MACHINE-to-Machine (M2M) communication or 
Machine-Type Communication (MTC) is expected to 
be one of the major drivers of cellular networks [1] and 
has become one of the focuses in 3GPP [2-4]. Behind 
the proliferation of M2M communication, the congestion 
problems in M2M communication become a big concern. 
The reason is that the device density of M2M communication 
is much higher than that in traditional Human-to-Human 
(II2H) communication [2,3]. For example, it is expected in 
[3] that 1,000 devices/km^ are deployed for environment 
monitoring and control. What is worse, in event-driven 
M2M applications, many devices may be triggered almost 
simultaneously and attempt to access the base station (BS) 
through the Random Access Channel (RACH) [5]. Such high 
burstiness can result in congestion and increase access delay, 
which motivates our research. 

Several efforts have been made in 3GPP to alleviate the 
radio congestion on RACH. Back to 3GPP LTE (Long Term 
Evolution), since the amount of exchanged data grows rapidly, 
congestion control has been addressed and access class barring 
(ACB) scheme has been proposed for overload protection 
[6]. In ACB scheme, devices are divided into several ac- 
cess classes. Before establishing a connection, the device is 
required to perform the ACB check and randomly transmit 



request packets with a probability broadcasted by the BS. 
Therefore, the traffic load can be reduced by choosing a small 
transmission probability. It is noticed that ACB scheme is a 
good candidate for congestion control in M2M applications, 
though modifications in accordance with the specific features 
of M2M applications is necessary [7]. Therefore, in [8], a two 
stage access control scheme, which consists of a primary level 
and a secondary level of access control barring, is introduced 
to provide prioritized M2M services based on their service 
attributes. Moreover, the authors in [9] propose a cooperative 
ACB scheme for balancing traffic load among BSs in a 
heterogeneous multi-tier cellular network. With cooperations 
among BSs, the congestion level can be reduced and the access 
delay can be significantly improved. However, a key problem 
arising in implementing these schemes is how to estimate 
the number of active devices and optimize the transmission 
probability. This problem becomes even worse in event-driven 
M2M applications which is characterized with highly bursty 
traffic. 

Essentially, the ACB scheme belongs to slotted-ALOHA 
(S-ALOHA) type schemes, which are widely applied for 
random access control. To address the instability issue of S- 
ALOHA [10], plenty of work has been done for deciding 
the protocol parameters and stabilizing the S-ALOHA system, 
which is briefly summarized in Section II. In these schemes, 
historical outcomes are applied to estimate the network status 
and adjust protocol parameters. Drift analysis is then used to 
design the schemes and prove their stability [11-13]. However, 
these schemes usually rely on the assumption that the traffic 
can be modeled as a Poisson process and only apply the 
observation in the previous slot for estimating the number of 
backlogged devices. Due to the burstiness, this assumption 
cannot be justified in the context of M2M applications and 
the observation in a single slot is not satisfactory for adjusting 
the protocol parameters in time. Thus, we try to make full 
use of the information provided by the historical outcomes 
for improving the performance of access control under bursty 
traffic. 

In this paper, we study adaptive S-ALOHA scheme for 
event-driven M2M communication and provide rigorous anal- 
ysis about its stability. As our main contribution, we propose 
a Fast Adaptive S-ALOHA (FASA) scheme for the random 
access control of M2M devices. A key characteristic of FASA 
is that the access results in the past slots, in particular, 
consecutive idles or collisions, are collected and applied to 
estimate the number of backlogged devices. This enables the 
fast update of transmission probability under highly bursty 
traffic and thus reduce the access delay. Furthermore, we prove 
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the stability of FAS A under bursty traffic. This is accompHshed 
by examining the T-slot drifts of the network status, which 
captures the memory property of FASA. Under interrupted 
Poisson traffic model [14], we show that the T-slot drifts of 
FASA have the required properties for stabilizing the scheme 
and the system is stable when the arrival rate is less than e^^. 
Numerical results demonstrate that using FASA scheme, the 
transmission probability of S-ALOHA can be adjusted in time 
and the access delay can be reduced to be very close to the 
theoretical lower bound under highly bursty traffic. 

The remainder of the paper is organized as follows. Sec- 
tion II summarizes the related work. In Section III, we present 
the system model, including the bursty traffic model for 
the event-driven M2M communications. In Section IV, after 
analyzing the limitations of traditional fixed step size policies, 
we propose the FASA scheme and design the parameters in the 
scheme based on drift analysis. Then in Section V, we study 
the T-slot drifts of FASA and prove its stability. In Section VI, 
simulation results are presented to evaluate the performance of 
the proposed scheme, compared with the theoretical optimal 
scheme in ideal case and two traditional adaptive schemes. 
Finally, we conclude our paper in Section VII. 

II. Related Work 

Radio Congestion Control in M2M communications: Due 
to the high equipment density, the radio congestion of M2M 
applications is still an open problem. Besides the efforts 
in 3GPP [7,9,15], there are a few publications addressing 
this issue. Since group-based feature appears in many M2M 
applications, some researchers propose hierarchical architec- 
tures, such as grouping scheme [16] and relay schemes [17, 
18], for alleviating the radio congestion on the RACH. In 
these architectures, group heads are selected for collecting 
the messages from the group members. However, efficient 
schemes for communications between the group head and 
members is required to reduce the access delay. At the same 
time. Adaptive Traffic Load Slotted Multiple Access Collision 
Avoidance (ATL S-MACA) mechanism proposed in [19] uses 
packet sensing and adaptive method to improve the access 
performance under high traffic load. However, the scheme is 
designed for Poisson traffic and is not suitable for event-driven 
M2M applications. 

Adaptive S-ALOHA Schemes: S-ALOHA is a fundamental 
scheme for ransom access control and the combinations with 
other techniques such as CDMA make it even more useful. 
However, the instability issue of S-ALOHA should be dealt 
with when being implemented in real networks. Two typi- 
cal classes of schemes, additive and multiplicative adaptive 
schemes, have been proposed for stabilizing S-ALOHA sys- 
tems [20]. In these schemes, the estimate of the network status 
is updated in a additive or multiplicative manner, respectively, 
and the transmission probability is adjusted accordingly. How- 
ever, as discussed in more detail later, traditional additive 
schemes such as Pseudo Bayesian ALOHA (PB-ALOHA) 
[21] estimate the number of backlogged devices based on 
the access result in the previous slot but cannot adjust the 
transmission probability in time under highly busty traffic, 
resulting in large access delay. On the other hand, because 
of the exponential increment in consecutive collision slots or 
decrement in consecutive idle slots, multiplicative schemes 
[22], e.g., Q-Algorithm in [23] and its enhanced version Q+- 
Algorithm in [24], can track the network status in a short 



period. However, the throughput suffers in these schemes due 
to the fluctuations in the estimation [22]. Therefore, we aim 
to design adaptive schemes that could adjust the protocol 
parameters fast under bursty traffic while retaining the same 
stable throughput as additive schemes. In our previous work 
[25], we propose a preliminary version of FASA based on 
some intuitive approximations and show its desirable prop- 
erties through numerical simulations. However, no rigorous 
analysis about the stability of FASA is presented in [25]. 

Drift Analysis for stabilization of S-ALOHA: Drift analysis 
is a theory for deducing the properties such as ergodicity 
of a sequence from its drift and is found useful in the 
design and analysis of adaptive S-ALOHA schemes [12,26, 
27]. The network status, which is represented by the number 
of backlogged devices and its estimate, could be viewed as a 
stochastic sequence. It is shown in [12] that when the drifts 
of the network status satisfies some criterions, the system is 
stable in the sense that the returning time can be bounded 
with high probability. Using the conclusion in [12], the most 
related work [13] studies the stability of PB-ALOHA scheme 
by defining a Lyapunov function to represent the network 
status and examining its drift. In all the work mentioned above, 
the schemes update the parameters based on the observation 
in the previous slot. Thus, the 1-slot drifts, i.e., drifts between 
two adjacent slots, are sufficient for studying the stability of 
the systems. When involving access results in multiple slots, 
however, the T-slot drifts are required to deal with the memory 
property of our scheme. Unlike 1-slot drifts, calculating T-slot 
drifts is non-trivial and we have to resort to approximations 
for obtaining their properties. 

III. System Model 

In this section, we describe random access control procedure 
for M2M communication as well as the traffic model, which 
will be used for studying the stability and evaluating the 
performance of the proposed scheme. 

A. S-ALOHA Based Random Access 

We consider a cellular network based M2M communication 
system for event detection. As shown in Fig. 1, the system 
consists of a BS and a large number of M2M devices. When an 
event is detected, certain number of devices are triggered and 
attempt to access the BS by sending request packets through 
a single RACH based on S-ALOHA scheme. 
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Fig. 1. S-ALOHA based access control of M2M communications 

The time is divided into time slots, each of which is long 
enough to transmit a request packet. Deferred first transmission 
(DFT) mode [21] is assumed, in which a device with a new 
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request packet immediately goes to backlogged state. In the 
tth slot, where t e Z+ := {0, 1,2,.. .}, all backlogged devices 
transmit packets with probability pt, which is broadcasted 
by the BS at the beginning of the slot. For the sake of 
simplicity, we assume that the request packets generated by 
M2M devices will eventually be transmitted successfully and 
a backlogged device will not generate any new requests since 
the new coming data can be transmitted as long as the device 
accesses the BS successfully. Moreover, an ideal collision 
channel is assumed, where the transmitted packet will be 
successfully received by the BS when no other packets are 
being transmitted in the same slot. 

The transmission probability pt is adjusted based on access 
results in the past. Let Zt denote the access result in slot t, 
and Zt = 0, 1, or c depending on whether zero, one, or more 
than one request packets are transmitted on the RACH. At the 
end of slot t, the BS decides the transmission probability for 
next slot based on the sequence {Z^, Zi, . . . , Zt}, i.e., 

Pt+i = ^tiZa,Zi, . . . , Zt). 

The objective of the BS is to maximize the throughput and 
minimize the access delay. It well known that, when Nt > 1, 
where Nt is the number of backlogged devices in slot t, using 
a transmission probability pt = l/Nt in slot t maximizes the 
throughput of the S-ALOHA system. However, the BS does 
not know Nt and has to obtain its estimate Nt based on the 
access results in the past. 

B. Traffic Model 

In order to capture the burstiness of event-driven M2M 
traffic, instead of traditional Poisson process, the arrival pro- 
cess is modeled as an interrupted Poisson process, which was 
suggested by Hayward of Bell Laboratories for simulating 
overflow traffic [14]. 

Interrupted Poisson process can be viewed as a Poisson 
process modulated by a random switch and will be discretized 
according to the slotted structure of the scheme. Let Yt and 
At respectively denote the number of events happening and 
the number of devices triggered in slot t. Assume that events 
happen independently and identically in each slot and at most 
one event happens in one slot, with 6 being the happening 
probability. Hence, in each slot t, PT{Yt = I) = 6 and 
PT{Yt = 0) = 1 — 6*. In addition, assume that the number 
of triggered devices follows Poisson distribution with mean 
A when an event happens, and no devices become active 
otherwise, i.e.. At ^ 'P{X) when Yt = 1 (ON-state), and 
At =i) when Yt = (OFF-state). Therefore, random variables 
{At : t £ Z+} are independently and identically distributed 
(i.i.d.) and the long term arrival rate can be calculated as 

A = Pr(yt = 0) •0 + Pr(yf = 1) • A = e'A. (1) 

In addition, the burstiness of the traffic is reflected by the 
variance of At, which is given by 

<t\ = 9{X + A^) - {eXf = A(l + A - A). (2) 

We note that the classic Poisson process is included as 
a special case of this model when 9 = 1. We will design 
and analyze adaptive S-ALOHA scheme based on this traffic 
model. Indeed, as will be discussed later, the scheme proposed 
in this paper could be stable under some other more general 
traffic models. 



IV. Fast Adaptive S-ALOHA 

The estimation of the number of backlogged devices plays 
an important part in stabilizing and optimizing the S-ALOHA 
system. In this section, using drift analysis, we first examine 
the limit of traditional fixed step size estimation schemes. 
Then, we propose and analyze a fast adaptive scheme, referred 
to as Fast Adaptive S-ALOHA. 

A. Drift Analysis of Fixed Step size Estimation Schemes 

Many schemes with fixed step size have been proposed in 
the literature to estimate the number of backlogged devices 
[20]. A unified framework of additive schemes is proposed 
and studied by Kelly in [11], where the estimate Nt is updated 
by the recursion 

Nt+i = max{l,iVt+ao/(^i = Q)+aiI{Zt = l)+aJ{Zt = c)}, 

(3) 

where qq, oi, and Qc are constants and I{A) is the indicator 
function of event A. 

With the estimation, the BS sets the transmission probability 
to = l/Nt for all backlogged devices, and thus the offered 
load p = NtPt = Nt/Nt, representing the average number 
of devices attempting to access the channel. To stabilize and 
optimize the S-ALOHA system, Nt needs to drift towards the 
actual number of backlogged devices Nt, especially when Nt 
is large. When Nt = n and Nt = n, the drift of the estimate 
can be calculated as follows [11]: 

E[Nt+i-Nt\Nt=n,Nt = h] 

= (ao - flc) 1 - - + (ai - flcj- 1 - - +ac 
\ n J n \ n J 

(ao - ac)e^'' + (ai - ac)pe^'' + Qc =^ 

as n — > oo, with n/h = p fixed. 

By properly choosing the parameters a; (i = 0, l,c) such 
that 4>{p) < if p < 1 and 4){p) > if p > 1, the 
estimate Nt will drift towards the true value and thus the S- 
ALOHA system can be stabilized. However, these fixed step 
size schemes are not suitable for systems with bursty traffic. 
When the estimate Nt deviates far away from the true value 
Nt, we have limp^o (pip) = and limp^oo 0(p) — o,c. These 
limits indicate that the drift tends to be a constant even when 
the deviation is large, which could result in a large tracking 
time. Thus, it is necessary to design fast estimation schemes 
for event-driven M2M communication. 

B. Framework of FASA 

As analyzed in the previous subsection, fixed step size esti- 
mation schemes such as PB-ALOHA may not be able to adapt 
in a timely manner for systems with bursty traffic because it 
always uses a constant step size even when the estimate is 
far away from the true value. We note that in addition to the 
access result in the previous slot, the access results in several 
consecutive slots will be helpful for improving the estimation 
as they could reveal additional information about the true 
value. Intuitively, collisions in several consecutive slots are 
likely caused by a significant underestimation, i.e., Nt ^ Nt, 
and the BS should aggressively increase its estimate. In 
contrast, several consecutive idle slots may indicate that the 
estimate Nt ^ Nt, and it should be reduced aggressively. 
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Motivated by this intuition, we propose a FASA scheme that 
updates Nt as follows: 

'max{l, Nt-l- ho{v){Ko,t A fc™)"}, if = 
Nt+i - <; Nu if = 1 

Nt + ^ + h,{y){K,^t^k,„T, if = c 

(4) 

where Xo.t and Xc.t are the numbers of consecutive idle and 
collision slots up to slot t, respectively; > 1 is an integer 
and K A k,n = min{K, fc„i}; > is the parameter that 
controls the adjusting speed; ho{iy) and hc{iy) are functions 
of u that guarantee the right direction of the estimation drift. 
In order to make the scheme implementable and its stability 
analysis tractable, we bound the update step size with k„i in 
this paper, which is different from that we proposed in [25]. 
However, the two schemes are almost the same as long as we 
choose a sufficiently large km- 

C. Design of hoiiy) and h(.{v) 

The functions hi^{v) and hc{v) are crucial for guaranteeing 
the convergence of the FASA scheme. Next we design hi^{v) 
and hc{v) by analyzing the drift of estimate Nf According 
to the structure of the proposed scheme, the evolution of Nt 
depends not only on the access result in the previous slot, 
but also results in the past km slots. Therefore, unlike the 
fixed step size schemes, accurate drift analysis is impractical 
for FASA because of its memory property. Thus, in order to 
make the problem tractable, we resort to approximation based 
on Lemma 1, which indicates the feasibility for approximating 
the distribution of access results in the past fc,„ slots with the 
network status in the previous slot. 

Lemma 1 For given e > 0, An, and An, there exists 
some M > 0, such that for any (n, h) G Hm, where 
Hm = {{n,n) : n > M ot n > M, n+An > 0, 71+A71 > 1}, 
we have 



(1- 



1 



\ n+An 



Af 



pe 



^I^+An-l 



An 



in 



(5) 

< e, (6) 



where p — n/n. 



Proof: See Appendix A. 
It is noticed that the distribution of the access result Zg in 
slot s is decided by Ng and Ng- In addition, for any given 
s e {t ~ km,t - fcm + 1, . . . , t - 1}, we have \Ns ~ Nt\ < 
km/{e - 2) + kl+''nmx{h„{i^),hc{iy)}, and \Ns - Nt\ is 
bounded with high probabiUty, i.e., PrdA^^ — Nt\ < B) ^ 1 
as B ^ 00. Thus, according to Lemma 1, when Nt and 
Nt are known and at least one of them is sufficiently large, 
the distribution of access results in the past km slots can be 
evaluated approximately, as well as the statistical characteris- 
tics of Kg^t and Kc.t- Therefore, in this section, we do not 
assume any knowledge about Ko t or Kc,t in slot t and will 
approximately calculate the drift of Nt conditioned on Nt and 
Nf In addition, we assume that Nt is large enough in the 
past km slots, and hence we can approximate max{l,a;} in 
(4) as X in the analysis later Rigorous analysis provided in 
Section V will show that the design with these approximations 
stabilizes the proposed scheme. 



Suppose that in slot t, the number of backlogged devices 
and its estimate are Nt ~ n, Nt — n, respectively, and thus 
the offered load p = n/n. When n or n is large, the drift of 
estimate Nt can be approximated as 



E[Nt+i^Nt\Nt^n,Nt^7\] 
^ ^ q,ip)E[ANt\{i,n,n)] 



(7) 



iG{0,l,c} 



where qaip) = e p, qi{p) = pe p, and qdp) = 1 - 
qo{p) — qi{p) are the probabilities of an idle, success, and 
collision slot, respectively; i?[AiVt|(i, n, n)] [i = 0,1, c) are 
the changes in Nt resulting from the corresponding updates. 

Obviously, E[ANt\{^in,h)] = since the estimated 
number remains unchanged when a packet is successfully 
transmitted in slot t. On the other hand, without memory 
about the access results in the past slots, K^^t and Kc,t 
are treated as random variables. Hence, E[ANt\{0,n,n)] and 
E[ANt\{c,n,h)] can be obtained based on the approximate 
distributions of A'o,t and Kc.t- 

First, to calculate the drift of estimate in an idle slot 
E[ANt\{0,n,h)], suppose that no packet is transmitted in 
slot t, then the estimated number will be reduced by 1 + 
ho{i'){Ko^t A kmY - Therefore, 

E[ANt\iO,n,h)] 
= - £ [l + /ioMfc^]Pr[A"o,t = fco|(0,n,n)] 



A;o = l 



[1 + /io(;^)C]Pr[A'o,t > km\iO,n,h)]. 



(8) 



Notice that /vo,t = fco(l < < km) holds when slots 
t — fco + 1> ^ ^ fco + 2, . . . , t — 1 are all idle while slot t — kg 
is not. Thus, for 1 < feg < fc,„, we have 

PT[Ko,t = ko\{0,n,h)] 
= PT[Zt-k„ 0\iZt-k„+i, - - - ,Zt,Nt,N^) = (0,...,0,n,n)] 

• n PAZs=0\{Zs+i,.-.,Zt,Nt,Ni)^{0,...,0,n,n) 

s—t — ko + l 

According to Lemma 1, when Nt = n or Nt = fi are 
sufficiently large, the distribution of access results in the past 
km slots can be approximated as that in the previous slot, i.e., 

VY[Zt-k,^Q\{Zt-k,+i,-.-,Zt,Nt,N't) = (0,...,0,n,n)] 
~ 1 - 90 (p), 

Vr[Zg^{)\{Zg+i,...,Zt,Nt,N't)^{{),...,Q,n,n)\ 
w qo{p), s = t - kf) + l,t - ko + 2, . . . ,t - 1. 

Consequently, for 1 < fco < km, 

PT[Ko,t = fco|(0,n,n)] « q^r'ip)[^ - loip)]- (9) 

Similarly, Ko^t > km holds if slots t— fc„i+l, fcm+1, 1 
are all idle and we can approximate the probability as 



Pr{Ko,t > km)~qt'\p) 



(10) 
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Substituting (9) and (10) into (8), we can calculate the drift 
of Nt in an idle slot as follows: 

E[ANt\{0,n,n)] ^ ~ J2 + ho{'y)k^oho"'\p){^ - M] 

ko = l 

-[l + /joMfc^]gS'""'(p) 
= -[1 + ho{v)fi{v,qo{p),k„i)], (11) 

where fi{iy,q,k„i) is defined as 

fi{l^,q,k^)^ J2 fcV"'(l-g) + (fcm)V'""\ (12) 
fe=l 

and qo{p), fcm) is the approximate expectation of (A'o,t A 
kmY conditioned on {Zt, Nt, Nt) = (0,n, n). 

Second, we can calculate the drift of estimate in a collision 
slot in a similar fashion as follows; 

E[ANt\ic,n,h)] 
« (e - 2)-i + h,{iy)E[{K\t A krny\{c,n, n)] 
= (e - 2)-i + K{v)ti{v, qaip), k^). (13) 

Therefore, the drift of estimate for FASA can be approx- 
imated by substituting the expressions of E[/S.Nt\{i,n,n)\ 
(i = 0, l,c) into (7): 

E[Nt+i~Nt\Nt = n,Nt=fi] 
~ -qo{p)[^ + ho{iy)p{i^,qo{p),k^)] 

+ qc{p)[{e - 2)-l + K{v)p{v, qdp), k^)] 

"^'vip)- (14) 

In order to keep the offered load p staying in the neighbor- 
hood of the optimal value p* = 1, it is reasonable to require 
that (^(1) = 0. In other words, letting q^ ~ 90(1) = e^^ and 
q* = qc{l) = 1 - 2e^\ we expect that 

(^(l) = -9o[l + hQ{v)p.{v,qQ,k^)] 

+ ql[{e - 2)-l + K{v)pl{v, q*,k^,)] 

= ~ho{iy)qQp{iy,q^,km) + hc{v)q*p{i',q*,k,n) = 0. 

(15) 

Hence, for given km > 1, to satisfy the condition in (15), 
we can select the following h(j{v) and hc{v): 

ho{v) ^ ■q[qap{y,ql,km)r^ , (16) 

hcif) =v[qcPi'^'(lc^l^-m)]~^, (17) 

where ?/ > is a constant. 

The chosen /lo(i^) and hciiy) guarantee that Lp{l) = 
and thus provide a necessary condition for FASA to track 
the number of backlogged devices. Furthermore, Theorem 1 
shows a desirable property of FASA, with which the estimated 
number Nt roughly drifts towards to the true value Nt, 
eventually yielding p = Nt/ Nt ~ I. 

Theorem 1 Given that hQ{v) and hc{v) are defined by (16) 
and (17), respectively, the approximate drift of FASA (p{p) is 
a strictly increasing function of p. In addition, ip{p) < when 
< p < 1 and ip{p) > when p > 1. 

Proof: See Appendix B. 
In order to understand better the behavior of the scheme, 
we now present the approximate drift of estimate for FASA 



with 77 = 1 and = 1, 2, 3. Assume that k„i is large, and thus 
the distribution of Ki t (* = 0, c) can be approximated as a 
geometrical distribution with success probability 1 — qi{p)- In 
addition, for ly S p{iy, qi{p), is approximately the i^th- 
moment of a geometrically distributed random variable with 
success probability 1 — qi{p), and its closed-form expression 
can be obtained. Consequently, the results obtained in our 
previous work [25] can be applied directly. Fig. 2 shows the 
approximate drift of estimation versus offered load, where 
the subscript of FASA represents value of (77, i/). It can be 
observed from the figure that when the estimated number 
deviates far away from the actual number of backlogged 
devices, i.e., p « or p ^ 1, FASA adjusts its step size 
accordingly, while PB-ALOHA still uses the same step size. 
Therefore, using FASA results in much shorter adjusting time 
than PB-ALOHA, and thus could improves the performance 
of M2M communication systems with bursty traffic. Note that 
the drifts of multiplicative schemes such as (5+-Algorithm are 
not illustrated here since they depend on not only the offered 
load p but also the estimate Nt- 




Fig. 2. Drift of estimation. 

V. Stability Analysis of FASA 

In this section, we use drift analysis to study the stability of 
the proposed FASA scheme. The M2M traffic is modeled as 
an interrupted Poisson process presented in Section III. In fact, 
as we will discuss later, the proposed scheme can be stable 
under other more general arrival processes. 

Unlike traditional adaptive schemes, the access results in the 
past consecutive slots are used in FASA to accelerate the speed 
of tracking, which makes it difficult to obtain the accurate 
drift of estimate. However, from the stability point of view, we 
concern mostly the scenarios where the number of backlogged 
devices or its estimate is large and hence approximation can 
be applied in these cases. To deal with the issues caused by 
the memory property of FASA, we analyze its T-slot drifts 
rather than the 1-slot drifts, which are introduced for ana- 
lyzing traditional adaptive ALOHA schemes [12,13,22,27]. 
By constructing a virtual sequence, we show that the T-slot 
drifts of FASA have the properties required for stabilizing the 
system as long as the number of backlogged or its estimate is 
sufficiently large, and these are similar to the properties of PB- 
ALOHA. Therefore, with slight modification, the Lyapunov 
function based method proposed for PB-ALOHA [13] can be 
used to prove the stability of FASA. 
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Consider the FASA scheme proposed in (4) under inter- 
rupted Poisson arrival process with average arrival rate A. We 
define a sequence Xt — {Nt, Nt, Kt), where Kt represents 
the memory of access results in the past consecutive slots and 
is defined as Kq = and for t > 0, 

-(A'o.t-i A fc,„), ifZt_i = 0, 

i^t = <! 0, if Zt-i = 1, 

Kc,t-i A km, if Zt-i = c. 

Recall that /^o,t-i and Kc^t-i is the number of consecutive 
idle and collision slots up to slot t — 1. Given initial value 
Xo = (0, 1, 0), each component of Xt evolves as follows when 

t>0 : 



denote the estimate error in slot t. Conditioned on the state of 
Xt, we define the T-slot drifts of Nt, Nt, and Nt as follows: 



Kt+i = < 



\Kt 



if Kt < 0, Zt 
if Kt > 0, Zt 
if Zt - 1, 
if Kt < 0, Zt 
if Kt > 0, Zt 



0, 
0, 

c, 
c, 



(18) 



(19) 



1) A kr, 

-1, 

0, 

1, 

[{Kt + l)Ak,n, 
Nt+i = max{0, Nt + At - I{Zt = 1)} 

rmax{l,iVt+i - 1 - h„{u)\Kt+i\''}, if Zt = 0, 
Nt+i = <^ Nt, if Zt = 1, 

{Nt + ^ + h,{v){Kt+iY, iiZt^c. 

(20) 

It is easy to verify that Xt = {Nt, Nt, Kt) is a Markov 
chain on a countable state space Sx- The main result of 
this section reveals the geometrical ergodicity [12] of Xt, 
which is described in Theorem 2. As pointed out in [12], 
the geometrical ergodicity is a weaker form of ergodicity and 
indicates the existence of steady distribution for each initial 
state. 

Theorem 2 If < X < e~^, km. > 1, haiv) and hdv) 
are given by (16) and (17), then the Markov Chain Xt is 
geometrically ergodic. 

Proof: The proof of Theorem 2 is based on the drift 
analysis. Specifically, the proof involves three steps, which 
are outlined as follows and presented afterwards: 

Step 1 - Approximation of drifts: To deal with the impact 
of the memory in the proposed scheme, rather than 1-slot 
drifts in the existing works, we study the T-slot drifts for 
our scheme, which are then approximated by constructing a 
virtual sequence X'tj^^ conditioned on the state of Xt in slot 
t. 

Step 2 - Property analysis of drifts: Based on the approxi- 
mation of T-slot drifts for FASA, we obtain the properties of 
the T-slot drifts required for guaranteeing the stability of the 
scheme. 

Step 3 - Stability analysis based on Lyapunov function: The 
Lyapunov function defined in [13] is adopted for the proposed 
scheme. Then with the the properties obtained in Step 2, we 
show the geometrical ergodicity of Xt by analyzing the drifts 
of the Lyapunov function. 

Step 1 - Approximation of Drifts 

In order to analyze the stability of FASA, we evaluate the 
change of Xt from slot t to slot t + T. Let Nt ^ Nt - Nt 



E[Nt+T 


~Nt 


Xt 


= {n 


h,k)] 


E[Nt+T 


-Nt 


Xt 


= (n 


n,k)] 


E[Nt+T 


-Nt 


Xt 




h, k)] 



dT{n, fi, k) — dxin, h, k). 



(23) 



It is difficult to calculate the drifts defined above and we try 
to obtain the properties of them by introducing an approximate 
version of Xt- Let {Z'tj^^ : s G Z+} be a ternary independently 
and identically distributed (i.i.d.) random sequence, whose 
distribution is given by 

P^iZ't+s = «) = « = 0, 1, c, 

where p = n/h. Then we construct a virtual sequence Xt^^ = 
{N^^g,N^^g,Kt_^_g) based on Xt and Z't^^ as follows: when 
s = 0, X[ ^ Xt = (n, h, k); when s > 0, K'^^^, N^^^, and 
Nt^s are updated in a similar way as (18) - (20), respectively, 
with Zt replaced with Zt_f_g. However, unlike the updates in 
Xt, we allow N^^^ to be negative and Nj-^^ to be less than 
1. 

Obviously, Xj.^^ is a Markov chain and its transition proba- 
bilities are fixed and determined by the state of Xt . We define 
its T-slot drifts d'rp{n, h, k), d'rp{n, h, k), and d'rp{n, h, k) sim- 
ilarly to (21) - (23). When T is given, we show in Lemma 2 
that the drifts of X^.^^ can be used to approximate the drifts 
of Xt from slot t to slot t + T, when either Nt or Nt is 
sufficiently large. 

Lemma 2 Given T > and e > 0, there exists some M > 0, 
such that if Nt ~ n > M or Nt — n > M, then 



\dT{n, h, k) — d'rp{n, n,k)\ < e, 
Idxin, n, k) — d'rj^{n, n,k)\ < e. 



(24) 
(25) 



Proof: See Appendix C. 
According to Lemma 2, for given T, the differences 
between the T-slot drifts of X't^^ and Xt+s can be made 
as close to zero as desired by letting n or n be sufficiently 
large. Next, we evaluate the drifts of X't^^ for obtaining the 
properties of drifts for FASA in Step 2. 

1) d'rp{n,n, k): Since in the virtual sequence Nj.^^ 
is allowed to be negative, we can easily have 



T-l 

d'r{n,h,k) = E[Y,{At+s - Z't+M - {n,n,k)] 

s=0 

= T{X-pe~n- 



2) d'rp[n,n, k): In slot t + s, the update of N't^g depends 
on both and Z't^^. Notice that the sequence K't_^_g is a 

Markov chain on a finite state space S/<- = {—km.,— km + 
1, . . . , 0, . . . , km}- Since the distribution of Z'tj^^ is fixed, by 
showing the ergodicity of K[_^g, we are able to approximate 
the T-slot drift by analyzing the stationary behavior of K'tj^^. 
Specifically, the transition of A'j^^, depends on the value of 
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I ^ and the 1-step transition probabilities is given by 



'qoif), if fc < 0,j = niax{fc - 1, -fc,„}, 

qo{p), if fc > 0, j = -1, 

qi{p), if i = 0, 

Qcip), iffc<0,i = l, 

qdp), if fc > 0,j = niin{fc + 

0, else. 



It is easy to verify that when p > 0, Kf_^_g is irreducible 
and aperiodic. Thus, Kf_^_g is ergodic and there is a unique 
stationary distribution. Now we study the stationary distribu- 
tion of A'j_|_g and the drift of Nl._^_g in the steady state. Define 
a 1 X {2km + 1) vector as follows: 

'^(p) = k-fe,„(p)>^-fc™+i(p)> • ■ ■,T^o{p), ■ ■ ■ ,7rfc„(p)], 



where the elements are given by 



7rfc(p) = < 



i'\p)[l 

nip), 



rc{p)[i 



- 90 (p)], 
(lc{p)], 



if fc = 

if -k,,^ 

if fc = 0, 

if 1 < fc < fc„„ 

TTj{p) = 1 and TTjip) 



1< fc < -1, 



1, 



iQTip)^ 
It can be verified that J2'i=-k 
J2k=-k^ ^k{p)Pkjip) for all fc e Sk- Hence, tt{p) is the sta- 
tionary distribution of K'f._^^. Using the expression of ir{p), we 
can verify that (p{p) defined in (14) represents the stationary 
drift of N^^s, which is the 1-slot drift of Nj-^^ when K'^^^ is 
in the steady state. Consequently, with the ergodicity of K[^g, 
we have 

1 . 

Vip)- 



lim —d'rr(n,h,k) 



Moreover, in order to use Lemma 2, we expect to find a 
common T such that (24) and (25) hold for some given e and 
for all (n, n, fc) e Sx, which requires the uniform conver- 
gence of ^d'rp{n,n, fc). In fact, by analyzing the evolution of 
i^t'+s, we can show that ^d'rp{n,n, k) converges uniformly 
in {n,n,k). First, by multiplying the transition probability 
matrix km times or analyzing the event that /Vj^j, = j, 
we can see that for any fc,j G Sk, the fc'm-step transition 
probability pl^j^^'if) = T^jif)- For example, for any fc e Sk 



and j e (-fc™,0), Kl 



t+k„ 



j holds if and only if Z^' 



t+s 







for all s = km — 1, hn — 2, fc,„ — j + 1, while Z^_^_,^^_- ^ 0, 

so ptj\p) = 'zl/k/o)!! - ^oIp)] = 7rj(p). Consequently, for 
any state fc G S^, we have 



Pr(-R't'+^ = j) = TTj{p), s > km, 

and thus when s > km the drift of A^t_|_g in each slot is exactly 
ip{p)- Then, with the fact that for any {n,n, fc) G Sx, 



d'k,^-i{nji, fc) - {km - ^)v{p) 
1 



< {km - 1) 



fc,''„max{/io(i'),/ic(j^)} 



we know that as T tends to infinity, ^d'j,{n,h,k) converges 
to (p{p) uniformly in {n,n,k). Thus, the difference between 
^d'rp{n, n, fc) and Lp{p) can be made as close to zero as desired 
by choosing a common T for all (n, h, fc) G Sx- 



3) d'rp{n,n,k): Since d'rp{Ti,h,k) ~ d'j^{n,h,k) — 
d'rp{n,n,k), we introduce the following function to approx- 
imate j;d'rp{n, h, fc): 

i;{p,\)^^{p)-{X-pe-n 

= + hc{>y)p{>^, 9c(p), fc„0](l - e-" - pe~P) 

e — 2 

- [1 + ho{v)p{v, qQ{p), km)]e^'' + pe^'' - A. 

With the uniform convergence of ^d'rp{n,n,k), we know 
that for any given A > 0, as T ^> oo, 

^d'rp{n,n,k) ilj{p,\) 

uniformly in (n, n, fc). 

Step 2 - Property Analysis of Drifts 

The evolution of estimate error Nt is critical for showing 
the stability of the scheme. We first show in Lemma 3 that the 
approximate drift '4'{P: A) has the same properties as those for 
the 1-slot drift of PB-ALOHA and then present the required 
properties of T-slot drifts of FASA in Lemma 4. 

Lemma 3 Given km as a positive integer, i}}{p,X) has the 
following properties: 

a) For any A, the function ip{p,X) is strictly increasing in 

P- - _ 

b) For any A G (0,e~^], there exists a unique p = w(A) G 

(0, 1], such that ip{p, A) = 0. 

c) IfXe (0,e-i), then w(A)e-"(^) > A. 

Proof: See Appendix D. _ 
Let /3 = a;(A) denote the root of tp{p, A) = for given A. 
Similarly to the method in [13], we partition the state space 
into the following four parts: 

S^^M = {{n,n,k) : n> M or 7i > M, 

TL 

n 

TL 

Km = {{n,n,k) : n > M, - < (3 - j,k e Sr}, 



Km 



{{n, n,k) : n> M, - > l+"f,k e Sk} 



^ {{n, n,k) -.rK M, n< M,k € §k} 



and let R^^m 

With Lemmas 2 and 3, we present the properties of the 
T-slot drifts in these regions in the following lemma. 

Lemma 4 There exist some 7 > 0, (5 > 0, T > 0, and M > 0, 
such that 57 < /3 and 

dT{n,h,k) < —TS, y{n,h,k) G ^5^^^/, (26) 

dT{n,h,k) < -TS, y{n,h,k) G K^m^ (27) 

dT{n, h, fc) > T(5, V(n, h, fc) G (28) 

Proof: See Appendix E. 
Intuitively, according to Lemma 4, when the estimate Nt 
is close enough to Nt, positive number of devices will access 
successfully and leave the network in the following T slots. 
On the other hand, the deviation of the estimate Nt from Nt is 
expected to decrease when it is larger than a certain threshold. 
These properties guarantee the stability of FASA, as presented 
in Step 3. 
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Step 3 - Stability Analysis Based on Lyapunov Function 

Lemma 4 shows that with a sufficiently large T, the T- 
slot drifts have similar properties to the drifts of PB-ALOHA. 
Hence, when observing the system every T slots, the Lyapunov 
function based method for PB-ALOHA can be used for 
analyzing the stability of FASA. Next, we provide an outline 
of using the Lyapunov function based method to prove the 
stability of FASA. For more details about this method, it is 
recommended to refer to [13]. 

Assume that T, M, 7, and 5 are fixed and that inequations 
(26) - (28) hold. We use the Lyapunov function defined in 
[13]: 



V{n, n, k) = max < n, 



1 + 37 
37 



/?-37 



(n — n) 



n, 

/3-37 
1-/3+37 



(n - n), 



" 1 - /J + 37 
if {n,n,k) £ 6*3^,1/, 
if {n,n,k) S R^j^m, 
if {n,h,k) e R^^^M- 



(29) 



We will show that if J is sufficiently large, there exists some 
A > such that 

E[V{Nt+jT,Nt+jT,Kt+jT) - V{NuNt,Kt) + A; 

\Tt\ < 0, (30) 

where J't is the a-field generated by {As-i, Ng, Ns, Kg : 
s < t} and for random variable X and event A, the notation 
E[X;A\T] stands for E[XI{A)\J^]. 
For given t >0 and integer J, let 

Tj = min{j > : ^ At+s > JT}. 

s=0 

Similarly to [13], we then analyze the drift of the Lyapunov 
function by considering the unlikely event {tj < J} and likely 
event {tj > J} separately. 

Using Chernoff bound [28], we can show that the following 
results also hold for interrupted Poisson process: 



lim Pr(T,/ < J) = 0, 

J— >oo 



lim E 

J— ^CXD 



JT 

hJT + h^A 

s=0 



t+s 



TJ<J 



(31) 



0, (32) 



where li, I2 are arbitrary given constants. Thus, for any 

(n, h, k) £ §A', as J — > cxD, we have 

E[\V{Nt+,jT,Nt+jT,Kt+,jT) ~ V{NuNuKt)\- 

TJ < J|Xt = (n,ri,fc)] ^0, (33) 

implying that this expectation can be made as close to as 
desired by choosing a sufficiently large J. 

Now consider the event tj > J. Based on the value of 
Xt = [n, h, k), we study the drift of the Lyapunov function 
in the following five cases: 
a) {n,h,k) e 5'2^.a/+(jt)2; 



b) (n, h, k) g i?+ 



■47, M+(JT)2 



c) {n,h,k) G Ri^^M+iJT)^' 

d) (n, n, k) G R4:j,m+{jt}^ H R2'y,M+{jTy 

f) {n,h,k) £ i?47,A/+(,7T)2 ^1 -R27,M+( JT)^ ■ 



In any of these cases, following the approach in [13], we 
can show that when J is sufficiently large, there exists some 
A > 0, such that inequation (30) holds. 

Take case a) as an example. According to Lemma 3.4 in 
[13], if Xt = {n,n,k) G 'S'27,a/+(jt)2, then we choose a 
sufficiently large J, such that {Nt+s, Nt+s, Kt+g) G S^^^m 
for all s = 0,1,..., JT, and 

V{Nt+jT,Nt+jT, Kt+jT) = Nt+jT, for all j G [0, J]. 

Thus, choosing a sufficiently large J such that Pr(r7 > J) > 
1/2, we have 

E[V{Nt+jT, Nt+jT, Kt+jT) - ViNt,Nt, Kt); 
Tj > J\Xt = {n,h,k)] 
= E[Nt+jT - Nt; TJ > J\Xt = [n, n, k)] 
.1-1 

= {n,n, k)] 

3=0 

STJ 



< ~5TJPt{tj > J)< 



(34) 



Combining (33) and (34), we know that there exists some J 
such that inequation (30) holds for some given A > 0. 

Now we are able to use the results about the hitting time 
bounds implied by drift analysis in [12]. Let 

1 + 37 /3 - 37 



A = [M + [JTy] max{l, ■ 



-}• 



37 'l-;3 + 37- 

Note that for any Kt G E>k, whenever Nt > M + {JT^ or 
Nt> M + {JTf, (30) holds. According to Theorem 2.3 in 
[12], for any initiate state, the returning time = min{t > 
: V{Nt, Nt, Kt) < A} is exponential type, which implies 
that Xt is geometrically ergodic and concludes the proof of 
Theorem 2. ■ 
Similarly to the discussion in [13], from the proof of 
Theorem 2, we know that the proposed FASA scheme is stable 
under more general traffics, as long as the average arrival rate 
A < and the traffic model satisfies the conditions in (31) 
and (32). 

VI. Simulation results 

In this section we evaluate the performance of the proposed 
scheme through simulation. We first examine the tracking 
performance and the effect of control parameters f and 77 on 
the access performance. We then study the access delay of the 
proposed scheme, including both the cases of single event and 
multiple events reporting. 

We compare the performance of our FASA scheme, the 
ideal policy with perfect knowledge of backlog, PB-ALOHA 
[21], and Q+-Algorithm [24]. With perfect knowledge of Nt, 
the ideal policy sets transmission probability at pt = 1/Nt 
for Nt > 0. Thus, the ideal policy achieves the minimum 
access delay of S -ALOHA and serves as a benchmark in the 
comparison. For PB-ALOHA, we use the estimated arrival 
rate At ~ e~^, as suggested in [13]. Q+ -Algorithm belongs to 
the class of multiplicative schemes which is first proposed by 
Hajek and van Loon [22]. In Q+ -Algorithm, Nt is updated as 
follows: 

Nt+i = max{l, {I{Zt = 0)/Co+/(^t = l)+Cc/(^t = c)]^*}, 
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where Co = 2°-25 « 1.1892 and Cc = 2' 
suggested in [24] for optimal performance 



0.35 



1.2746 are 



A. Performance of tracking 

In order to gain more insights into the operation of the 
estimate schemes, we treat adaptive S-ALOHA schemes as 
dynamic systems and study their step responses, where the 
number of backlogged devices Nt =Q when t < and Nt = n 
for all t > 0. We examine the tracking performance of schemes 
for n = 500, 1000, and 2000. For the sake of simplicity, we 
fix the value of km in FASA at 20, since the effect of km 
vanishes when it is large enough due to the the exponential 
decay of distribution of Kt ^- 

Before quantitive analysis, we first show the evolution of 
estimations under different conditions, which gives us some 
perceptual understanding about the behavior of these schemes. 

Fig. 3 shows the evolution of Kt and Nt for different 
adaptive schemes, where the subscripts of FASA represents 
values of rj and v. Average values of Kt and Nt in each 
slot are obtained from 4000 independent experiments. It can 
be seen from this figure that, unlike the almost linearly 
increasing in PB-ALOHA, the estimate Nt given by FASA 
increases slowly at the beginning, but speeds up due to the 
consecutive collisions, i.e., the increment of Kf. When the 
estimate gets close to the true value, success and idle slots 
occur and hence the increment of the estimate slows down. The 
estimate of Q+ -Algorithm follows the same trend as FASA 
and speeds up even faster on average than FASA because of 
the exponentially increment. Though the drift of the estimate 
has the right direction, Q+ -Algorithm turns out to be a bias 
estimate scheme, which will result in the suffering of the 
throughput at steady state. When comparing the curves of 
estimate for FASA with different parameters, we can see that 
with larger rj or i/, the estimate adjusts faster. 

Fig. 4 shows the evolution of Kt and Nt for different 
numbers of devices n. As shown in the figure, for larger n, 
there are more consecutive collisions at the beginning, which 
results in larger increment of estimate. After the peak point, 
the average value of Kt mainly depends on the offered load 
Pt = n/Nt or the ratio of Nt/n, so does the step size. For 
example, since iVio ~ 300 for n = 500, and N^o w 1200 
for n = 2000, i.e., the ratios of Nt/ri in these slots for 
n = 500 and 1000 are both about 0.6, they have almost 
the same average value of Kt, which is about 1. Though Kt 
decreases and FASA behaves like a fixed step size scheme as 
the estimate gets close to the true value, the value of Kt with 
larger n is larger than that with small n. Hence, we can expect 
that the time taken by FASA to catch up the true value or 
certain proportion of the true value will increase more slowly 
than linearly in n. 

Because the access delay depends on the throughput, two 
throughput-oriented metrics are introduced to measure the re- 
sponse speed and stationary performance: Q%-x% throughput 
rising time and stationary throughput. The 0%-x% throughput 
rising time is defined as the time required for the expected 
throughput to rise from 0% to x% of the optimal value e^^. 
For X = 10, 50, and 90, they are equal to the time required 
for the estimated number of backlogged devices Nt to rise 
from 0% to 20.45%, 37.34%, and 65.25% of the true value n, 
respectively. Stationary throughput is the average throughput 
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Fig. 3. Evolution of estimation with different parameters (n = 1000). 



after the time that the expected throughput reaches 90% of the 
optimal value. 

As shown in Table I, for the same x, the Q%-x% rising 
time (unit: slot) of PB-ALOHA almost linearly increases in 
n and is much larger than that of Q+-Algorithm and FASA. 
For instance, when n = 1000, the 0%-50% rising time of 
FASA with r/ = 1 and = 2 is about 1/20 of that of PB- 
ALOHA. Moreover, due to the aggressive update in FASA, 
the Q%-x% rising time increases more slowly rather than 
linearly as the number of devices n increases, especially when 
X is small. Comparing the rising time of FASA with different 
r] and ly, we observe that the increment of 7] or ly results 
in reduction of rising time. With multiplicative adjustment, 
it takes longer time than FASA for Q+-Algorithm to reach 
10% of the optimal value but shorter time to increase the 
expected throughput from 10% to 90% of the optimal value. 
In addition, the increment of the rising time in Q+ -Algorithm 
is tiny as n grows. However, aggressive adjustment of estimate 
usually results in large fluctuation at the steady state, and thus 
lower stationary throughput, which is shown in Table II. Thus, 
trade-off between the rising time and stationary throughput is 
necessary for choosing the values of parameters. 

B. Access delay 

In this section, simulation results about the access delay 
of adaptive S-ALOHA schemes are presented. In some event- 
driven M2M applications, response can be taken with partial 
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Fig. 4. Evolution of FASA for different n (r] = I, u = 2). 



TABLE I 

0% - X% THROUGHPUT RISING TIME 



n 


X 


PB-ALOHA Q+-AIg. FASAi. 


1 FASAi 3 


FAS A 1,2 


FASAi/2,2 


FASA2 2 




10 


59.8 


21.0 


9.1 


6.0 


7.1 


8.1 


5.0 


500 


50 


117.9 


23.5 


14.1 


8.0 


9.1 


11.8 


7.5 




90 


313.1 


27.5 


43.9 


14.5 


21.3 


31.9 


15.6 




10 


118.5 


23.0 


13.2 


7.1 


8.3 


10.1 


7.1 


1000 


50 


237.0 


26.6 


21.6 


9.4 


12.0 


15.0 


9.0 




90 


627.4 


30.7 


83.2 


24.8 


38.1 


60.0 


20.6 




10 


236.0 


26.0 


18.3 


8.1 


10.1 


12.1 


8.1 


2000 


50 


473.8 


29.5 


33.7 


10.3 


14.9 


19.7 


11.9 




90 


1254.5 


33.9 


165.1 


33.2 


62.5 


115.1 


38.1 










TABLE II 














Stationary throughput 






n 


PB-ALOHA 


Q+-Alg. 


FASAi,i 


FASAi.a 


FASAi _ 2 


FASAi/2 2 


FASA2,2 


500 


0.3686 


0.3529 


0.3679 


0.3656 


0.3670 


0.3676 


0.3653 


1000 


0.3684 


0.3521 


0.3678 


0.3663 


0.3675 


0.3681 


0.3668 


2000 


0.3681 


0.3523 


0.3678 


0.3674 


0.3677 


0.3678 


0.3675 



messages from the detecting devices and not all devices need 
to report an event. Thus, both the distribution of access 
delay for single event reporting and the long term average 
delay for repetitive event reporting are evaluated to study the 
performance of the proposed scheme. 

7 ) Single event reporting: We focus on the scenarios where 
a large amount of devices are activated to report a single event 
and study the distribution of access delay of different adaptive 
schemes. Assume that A^o = n devices are triggered at the 



same time when an event is detected, and attempt to access 
the BS on the RACH. The scenarios with the number of active 
devices n = 500, 1000, and 2000 are studied. 

Table III provides the y% access delay (unit: slot), which 
is the access delay achieved by y% of the active devices, and 
y is set to 10, 50, and 90. From Table III, we can see that 
the performance of the proposed FASA scheme is close to 
the benchmark with perfect information. For PB-ALOHA, it 
takes a long time to track the number of backlogged devices 
and few devices can access successfully during this period. For 
instance, the 10% delay is about two times of that for other 
schemes. For example, when n = 1000, the 10% delay of 
FASAi,2 is 290.8 slots while it is 542.4 slots for PB-ALOHA. 
With multiplicative increment, the Q+ -Algorithm can track the 
number of backlogs in a short time because of the exponential 
increment due to the consecutive collision slots. However, it 
takes longer for all the devices to access the channel under 
Q+ -Algorithm than under FASA due to the large estimation 
fluctuations in Q+ -Algorithm. Comparing the access delay 
achieved by FASA with different 1] and v, we observe that 
the 10% access delay is slightly smaller for larger r] or v, 
since they provide a quicker response ability. However, the 
larger fluctuation makes the 50% and 90% access delay for 
larger 77 and v close to, or even larger than that for smaller r] 
and V. 



TABLE ni 

j/% Access Delay 



n ?/ Peif.Info. PB-ALOHA Q+-Alg. 


FASAi, 1 


FASAi,3 


FASAi,2 


FASAi/2 2 


FASA2,2 


10 136.1 


271.9 


151.2 


155.2 


146.2 


146.4 


152.9" 


143.1 


500 50 680.6 


822.3 


712.5 


702.7 


694.4 


691.3 


700.6 


691.0 


90 1223.1 


1433.0 


1282.9 


1250.7 


1252.1 


1244.0 


1254.2 


1247.8 


10 267.3 


542.4 


298.6 


306.4 


285.9 


290.8 


303.3 


282.9 


1000 50 1351.9 


1648.2 


1426.0 


1394.3 


1384.7 


1385.7 


1385.8 


1378.5 


90 2440.8 


2871.4 


2558.7 


2496.0 


2488.2 


2484.1 


2489.9 


2483.5 


10 545.3 


1083.8 


585.1 


614.5 


563.7 


568.6 


592.4 


564.3 


2000 50 2712.1 


3294.9 


2855.8 


2793.4 


2749.2 


2752.1 


2779.1 


2744.6 


90 4886.3 


5745.8 


5124.3 


4983.1 


4944.9 


4944.0 


4968.6 


4934.5 



2) Repetitive events reporting with interrupted Poisson traf- 
fic: The events happen sequentially in the real system and we 
now study the long term average delay of adaptive schemes 
under interrupted Poisson traffic with different arrival rates and 
bursty level. In addition to average delay, we also define the 
normalized divergence as follows to quantify the divergence 
from the theoretical optimum performance: 



where D is the average delay of a particular scheme and D* 
is the theoretical optimal delay with perfect information. As 
pointed in the single event reporting case, the performance of 
FASA with different 77 and v are rather close. Thus, only the 
performance of FASA with ?7 = 1 and v = 2 is presented 
here. 

Fig. 5 compares the average delay and the normalized 
divergence of adaptive schemes under different arrival rates 
and fixed ON-probabiHty 9 = 0.0001. From Fig. 5(a), we 
observe that both PB-ALOHA and FASA are stable when 
the average arrival rate A < e^^ and experience finite access 
delays. For Q+-Algorithm, however, when the arrival rate is 
larger than about 0.352, the access delay grows unbounded, 
indicating that the algorithm with the given parameters is 
unstable for some A < e^^. As pointed out in [12], the 
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parameters in Q+-Algorithm should be carefully chosen to 
stabilize the scheme according to the value of A, which is 
not required in either PB-ALOHA or FASA. As shown in 
Fig. 5(b), when the arrival rate is close to zero, the divergence 
of Q+-Algorithm and FASA are larger than that of PB- 
ALOHA due to the fluctuation of estimation, while all the 
delays are very small. As the average arrival rate increases, the 
divergence of FASA decreases and gets close to the optimal 
value, since the estimate error becomes relatively smaller 
compared to the increasing number of backlogged devices in 
the system. For A larger than 0.1, The divergence of FASA is 
about 2.5%, while it is about 22% for PB-ALOHA. 

We point out that since the average access delay could be 
dominated by the time waiting in the system after the estimate 
catches up the true value, the improvement of performance 
by FASA does not seem to be significant from the average 
delay point of view. However, as discussed in the single 
event reporting cases, the 10% access delay can be improved 
significantly by FASA, which is very important to the event- 
driven M2M communications. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
Average arrival rale 

(a) Average delay (9 = 0.0001) 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
Average arrival rate 



(b) Divergence (6* = 0.0001) 
Fig. 5. Average access delay with different arrival rates. 

In order to examine the impact of burstiness, we present 
in Fig. 6 the divergence of access delay versus variance of 
the arrival process cr^ for A = 0.05 and 0.35. For the light 
traffic scenarios with A = 0.05, when the bursty level is low, 
i.e., the variance is small, the access delay obtained from all 
these scheme are close to the optimal value. The reason is 
that there is usually only one device is triggered in one slot 



when the estimate is usually set to 1 after several idle slots. As 
the traffic becomes more bursty, the divergence first increases 
owed to the rising time and fluctuation of estimate; and then 
the divergence of FASA and Q+ -Algorithm decreases for high 
bursty traffic because with aggressive update, they are able to 
track the status of the network quickly while the fluctuation 
become relatively smaller compared to the total number of 
backlogs. For the high traffic scenarios with A = 0.35, the 
divergences keep decreasing as the bursty level increase, while 
the proposed FASA scheme performs better than both PB- 
ALOHA and Q+-Algorithm under high bursty traffic. Since 
the traffic in the event-driven M2M applications is bursty, we 
believe that our proposed scheme will perform well for these 
applications. 



70 




10"' 10° 10' 10^ 

Variance 

(a) A = 0.05 





Variance 



(b) A = 0.35 

Fig. 6. Divergence of average access delay under different burstiness. 



vn. Conclusions 

In this paper, we proposed a FASA scheme for event- 
driven M2M communications. By adjusting the estimate of 
the backlogs with statistics of consecutive idle and collision 
slots, a BS can track the number of backlogged devices more 
quickly. That is a main advantage compared to fixed step size 
additive schemes, e.g., PB-ALOHA. Moreover, we studied the 
stability of the proposed FASA under bursty traffic, which is 
modeled as an interrupted Poisson process. By analyzing the 
T-slot drifts of the FASA, we showed that without modifying 
the values of parameters, the proposed FASA scheme is stable 
for any average arrival rate less than e~^, in the sense that 
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the system is geometrically ergodic. This property results in 
a much better long term average performance under heavy 
traffic loads, as compared with that of multiplicative schemes. 
In summary, the proposed scheme is an effective and stable S- 
ALOHA scheme and is suitable for the random access control 
of event-driven M2M communications as well as other systems 
characterized by bursty traffic. 

Appendix A 
Proof of Lemma 1 

Recall that for An ~ An = 0, it has been proved in 
[22] that when either the number of backlogged devices n 
or its estimate n is sufficiently large, the idle and success 
probabilities (and hence as well as the collision probability) 
can be approximated as functions of the offered load p ~ n/n. 
We generalize the results to any given Ati and An by showing 
that as n or n tends to infinity, the difference between the 
offered loads (n + An)/(n + An) and p can be ignored and 
the distribution of access results can still be approximated by 
the same functions of p. 

First, consider inequation (5), which is used for approxi- 
mating the probability of Zt = 0. Notice that 



Consequently, combining all the cases analyzed above and 
letting = max{M^°\ m[°I , p^°^ m[°^} follows that 

Next, we turn to the proof of inequation (6), which is 
used for approximating probability of Zt = 1. Similarly to 
inequation (5), 



< 



n 

—e 

n 

n 

— e 

n 



n - 
n - 



A?7. 

An 



(1- 



n 



An 



n+An-l 



An 



An 



An 



An 



(1 



\ri+An-l 



Using the result in [22], we know that there exists some 



M. 



> 0, such that G^^ < e/2 for any (n,n) G H 



M. 



(1) 



(1) 



< e/2 under certain 



(1 



1 



\n+An 



An' 



< 



(1 



hAri 



An' 



For according to Proposition 2.1 in [22], we know 



that there exists some M2^^ > 0, such that G'^' < e/2 for any 
(n, n) e H^jiQ), where Hj^^m = {{n,h) : n > M^^ or n > 

An > 0, n + An > 1}. Hence, we just need to 



(0) 



and we only need to show that 
conditions. 

Since lim„_>o pe~P = limp_>oo pe^'' = 0, there exist some 
p[^^ and such that < p[^^ < 1< and pe'P < e/4 
for any p satisfying < p < p\^^ or p > pj^'. Given p^^' £ 
(0, pj^'') and pj^-* G {p'^\oo), we study the following three 
cases based on the range of n/n: 

1-a) < n/n < p['\- Similarly to the analysis of O-a), 
since p']^^ < p\^\ there exists some m[^1 > 0, such that if 
n > m[^1, then (n + An)/{n + An) < p"^^ and hence 



(0) 



show that G^°' < e/2 when either n or n is sufficiently large 



G« < 

n 



n - 



An 



An 



< e/2. 



1-b) n/n > p^\- Similarly to cases O-a and l-a), we can 



Since limo 



there exists some p'"' > 1 such show that there exists some m[^2 > 0' such that if n > M| 



that e^P < e/4 for any p > />'•"> 
we analyze the value of G^^"* 
n/h < p^^\ respectively. 



f(i) 



. Given a number > p^^\ 
when n/h > p'^"^ and < 



then g[^\ < e/2 

have 

gW ^ ?e-"/ 



1-c) pj^-* < n/n < p'^\- When pj^'' < n/n < 



O-a) n/h > p^^\- When n/h > p^*^-* > p(°\ it is easy to 
show that there exists some m[^1 > 0, such that if n > m[^1, 
then (n + An)/(n + An) > (n - |An|)/(n + |An|) > p^°\ 
Hence, when n/h > p'°-' and n > m[^1, we have 

G[°^ < + e~^^ < e/4 + e/4 = e/2. 
0-b) < n/n < p^^\- Since n/n > 0, we have 



1 - 



_^ n n 
n h - 
1 + An/n 
1 + An/n 



At 



Ah 



cxp 



exp 



lAi 



nAh — hAn 
h{h + Ah) 
— hAn 



h{h + An) 



Since 



Gi' 



(0) 



1 — cxp(- 



n 



An, 



< 



1 — exp 



n 

nAn - 



n + An' 
nAn 



1 + An/n 



An/n 



- 1 



|An|/pW + |A 



< 



Ar 



0, 



CO, we have lim^. 



An) 



When < n/n < p("), 
nAn — nAn 



< 



n(n 
we have 
p(")|An| + |An| 



as n 

linift^oo i + An/fi *^^P["ri(rl+A^)"] 



l+An/n _ 



1 and 



(n + An) 



as 71 



n(n + An) 

oo. Therefore, limft^txj exp[ 

r(0) 



] = 1, and 

hence, there exists some M\^2^ such that Gi'' < e/2 for any 



nAh— nAn 
n(n+An) 
(0) 



(n, n) satisfying < n/n < p'"' and n > M} 



(0) 



•oo l+An/n 

1. Hence, there exists 

a m}^3 such that G^^'' < e/2 for any (n, n) satisfying 

pf ^ < n/n < p^^^ and n > Affg^ 

Therefore, G^^) < G^^' +G^^' < e for any (n,n) G i?M(i)' 
where 

M(i) = max{M(^\MW,Ml^\\p«AfW}. 

Finally, from the above analysis, choosing M = 
max{M(°\ M(i)}, we know that (5) and (6) hold for any 
(n, n) G Hm and this concludes the proof of Lemma 1. 
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Appendix B 
Proof of Theorem 1 

The proposition can be proved by calculating the derivative 

of Lp{p). 

For given values of v and fc,„, let 

JO) 



v'''Kp) = Mp)mNtmn,h)]\ 

= qo{p)[^ + ho{i')p{i',qoip),km)], 



and 



ip^'\p)^\qc{p)E[ANt\{c,n,h)]\ 

= Qc{p)[{e - 2)"^ + hc{v)p{v, qc{p),k^)]. 

Then = (p'^'^Hl) = + and ip{l) = -(^(o)(l) + 

(^(^) (1) = 0. Next, we claim that, for given > 0, q, fc,„) 
defined in (12) is an increasing function of q (0 < q < 1). 
This is because 

a A;,„-l 

? = E k^q^-'m -q)-l] + {k^ - l)k 
^ k=l 
k' 



J2k-q'-'[kil-q)-l] 



fc„-l 




>ik*) 



q)-l] + (/c™ - Ij/c^q'^ 



fc(l - g) - 1] + (fc™ - 



A:=l 



0, 



where A;* = [(1 — g) is the largest integer not greater 
than (1 - q)-\ and thus k" < {k*Y if 1 < fc < fc* and 
k" > {k*Y if fc > fc*. In addition, the idle probability 
9o(p) — e"'' is nonnegative and strictly decreasing in p. 
Hence, p{h',qo{p),km) is strictly decreasing in p and ip^^\p) 
is a strictly decreasing function of p. On the other hand, since 
Qcip) — 1 — e^'' — pe^P is nonnegative and strictly increasing 
in p, we can similarly show that ip^'^^{p) is a strictly increasing 
function of p. Thus, ip{p) = —(p'-^^p) + (p'-'^^p) is a strictly 
increasing function of p. Consequently, Lp{p) < ip{l) = 
when < p < 1 and Lp{p) > ip{l) = when p > I. 

Appendix C 
Proof of Lemma 2 

Because of the similarity, we present a complete analysis 
of inequation (24), while discuss briefly about inequation 
(25) at the end. Recall that we assume the same realization of 
arrival process At+s for Xt+s and Xl._^_^. In addition, given 
T and {n, n, fc), the number of departures between slot t and 
t + T is bounded and the number of new arrivals can also 
be bounded with high probability. In order to use Lemma 
1, we consider separately the events Y^^^o ^t+s > Ba and 



E 



T-l 



< Ba. 



1) E^Zo ^*+.^ > Ba 



Using Chernoff bound, we can show that the probability 
Pr(^^Jg At+s > Ba) decays exponentially as Ba grows 
to infinity. With the assumption of same reaUzation of At+s, 



the differences between the T-slot drift of Nt+s and Nt_^,^ is 
bounded by T, i.e., IdTirijh, k) — d'rp{n,n, k)\ < T. Therefore, 
for any given e > 0, there exists some Ba > such that 



E 

s=0 



At+s > Ba) < e/2. (36) 



Idrin, n, fc) — d'rp(n, fi, fc)|Pr( 
2J EL"o' At+s < Ba 

Let At^T = {At,At+i,...,At+T-i), Zt,T 
{Zt, Zt+i, . . . , Zt+T-i), and Ti'^j. — {Z^, Z[j^^, . . . , Z[^rj,_^). 
The set of possible values of {At^T,'Zt,T) and {At^T,'Z't t) 
are the same, denoted by Q,. Since Zt+s or Z't^^ has three 
possible values and YI^Zq -^t+s < Ba, is a finite set. 
Each pair (a, z) = (ag . . . a^-i, zo ■ • ■ ^t-i) G ^ results 
in corresponding drifts in both Nt+s and N't^^, denoted by 
AiV(a,z,„,ii,fc) and AiV('^^^^_-_^), respectively. 

One of the differences between Xt+s and Xj^, is that we 
do not limit the values of N'^^^ and N't+g, i-C-, ^t+s < and 
N't^g < 1 are allowed in the virtual sequence, which is not the 
case in Xt+s- However, noticing the fact that when n > T, 
the drifts of Nt+s and N|._^_g are only decided by (a, z), we 

for any (a, z) G O in this 



have AN, 



(a,z,n,n,A;) 



(a,z,n,n,/c) 



case. We first study the case where n >T, and analyze later 
the other cases where n is not large enough. 

Conditioned on 'Y^i^o ^t+s < Ba, we define the following 
probabilities: 



/A(a)=Pr(At,T = a), 

/(A,z)|x(a, z|n, n, fc) 
= Pr[(At,T, Zt.r) - (a, z)\Xt = (n, h, fc)], 

f{A,Z')\X'{a,z\n,n,k) 
= Pr[(Ai,T, Kt) = (a, z)|Xt' = {n, n, fc)], 

/z|A,x(z|a,n, n, fc) 
= Pr[Zt^T = ■z\M,T = a,Xt^ {n,n,k)], 

/z'|A.X'(z|a, ?^,"-,fc) 
= Pr[Zj,j, = z|At,T = &,X't = {n,n,k)]. 
Thus, when n >T, we have 



T-l 



\dT{n, fi, fc) — d'rp[n, h, fc)|Pr( 



s=0 



A 



t+s 



<Ba) 



< 



E [/(A,Z)lx(a,z|n,n,fc) 

(a,z)eO 

- ./(A,z']A'' (a, z\n, n, fc)] AA^(a,z,„,ft,fc) | 



(a,z)eO 



[/z|A,x(z|a,n,n,fc) 



- ./z'|A,A-(z|a, n, n, fc)]/A(a)AA^(a,z,n,ri,fc) I 

< C" E |/z|A,A'(z|a,ri, fi, fc) - /z'|A,x(z|a,n, fi, fc)| , 

(a,z)Gf2 

where C is the maximum value of |/A(a)AA^(a z, n.fi.fc) I for all 

1 J kjji } . 



fc^i 



1, 



(a, z) £ 51 and fc G {— fc 
Letting 

5z,(2s|a, zo ■ • • Zs-i) = Pr[Zt+s = 
Af^T ~ a,Zt . . . Zt+s-i 
.gz^(zs|a, Zt... zt+s-i) = Pr[.^i'+s 
Af,T ==a,Z;...Zj'+^_i 



..,0,. 

Zs\Xt 



k 



{n,n,k), 

Zq... Zs-l], 

Zs\X't = {n,n, fc), 

Zq... Zs-l], 
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we have 

/z|A,A'(z|a,n, n, fc) 

T-l 

= Pr(Zi = zq) gz,{zs\a, zq . . . z^-i) 

s=l 

T-l 

= ]^ Pr[^t+s = 0s|(7Vt+s,-/Vt+s) = (ns,ns)a,z,n,n,fc], 
s=0 

/z'|A,x(z|a,n,n,fc) 

T-l 

= Pr(^t = 2o) Jl gz'^{zs\a,zo . . . Zs-i) 



T-l 



n p^[^*' 



t+s 



a.z,n,n,j 



s=0 



where {ns,fis)a.z.,n,h,k is the number of backlogged devices 
and its estimate in slot t + s, given the initial state Xt = ~ 
(n, n, fc), At,T ~ a, and Zt.x = Zj = z. 
According to the construction of X[^g, the distribution of 
is fixed, i.e., Pr(Zf_|_j, = i) = qi{p)\ while for Xt+s, 
there are finite possible values of Nt+s ~ n and 7Vf+s — n, 
for all s = 1,2,..., T - 1, since the initial value of Kt is 
from a finite set. Therefore, using Lemma 1, we know that 
there exists some sufficiently large M[, such that if n > M[ 
OT n> M[, then for all s = 1, 2, . . . , T - 1, and (a, z) e Q., 



Vl[Zt+s = Zs\{Nt+s,Nt+s) = (?^s,'^s)a,z,n,r^,fc] 

where |0| is the number of elements in $7. Hence 

|/z|A,x(z|a,n,n,/i:) - fz'\A,x{A^,n,n,k)\ < 



< 



2TC\n\ 



2C\n\ 



<BA)<e/2. (37) 



and thus, 

T-l 

Idrin, n, k) — d'rp{n, n, fc)|Pr(^^ At+s 

s=0 

Therefore, combining (37) with (36) implies that (24) holds 
when n >T, and either n > M[ ot n > M[. 

Now consider the cases where n < T. In these cases, 
(At TiZt x) = (a, z) is an impossible event for some 
(a, z) G il, if it results in some Nt+s < 0. For these 
(a,z), we set AiV(a,z,„,ri,fc) = '^Nls,,zji.n.kr which does 
not affect the calculation of drifts since the probabilities of 
these events in Xt+s are zero. Notice that for these (a, z), 
there is at least one component of z equal to 1. Because 



lim 



0, we can find a number pi S (0, 1), such that 



pe-P < e/{2TC\n\) for all p = n/n < pi. Then the differ- 
ence between /z|A.x(z|a, n, n, k) and /z'|A,x(z|a, n, n, k) is 
bounded by e/{2TC\^\) and the same conclusion holds when 
n <T and n/fi < pi. 

Consequently, from the analysis above, we know that in- 
equation (24) holds when n > Mi or h > Mi, where 
AIi^max{T/pi,Mi/pi}. 

Finally, we discuss about inequation (25). For given 
(a, z) e ft, denote the corresponding T-slot drifts of Nt+s 
and Nt+, by A7V(a,z,„,fi,fc) and AiV^'^ ,^ respectively. 
Note that for inequation (24), we do not specially treat the 



case where n is not large enough (but n is large enough), 
since we still have AN(„ , „ i.\ = AN', . , x and the 

^a,z,7i,n,/i.; (a,z,7i,n,/e) 

update of FAS A guarantees that Nt+s > 1- However, 

AA^(a,z,n,«,fc) 7^ ^^(a,z,n,«,fc) ^OmC (a, z) 

when h is small. Since for aU (a, z) G il, both AiV(a ^ n^^) 
and ATV^'^ z n n fc) bounded uniformly in {n, n, k), the im- 
pact of AiV(a „ ,-,^fc) 7^ A7V('^^^^^ ^.^ can be made ignorable 
by making the probability of this event as close to zero as 



possible with the fact that limp_ 



0. Hence, similarly 



to inequation (24), we can analyze the following three cases 
to obtain the threshold for inequation (25): 

a) n>T,n>T[l + k';^ho{iy)] + l; 
h) n <T,h>T[l + C/io(i^)] + 1; 
c) n>T,n <T[1 + k-^^hoiiv)] + 1. 

Therefore, the proof of this lemma can be concluded by 
choosing M as the lager threshold for inequation (24) and 
inequation (25). 

Appendix D 
Proof of Lemma 3 



^(1 



a) Let 

"^'^Hp^A) = -er''ho{iy)p{iy,qo{p),k„ 



pe- 



e-P - pe'P > 0, 



+ (l-e P - pe P)hc{v)p.{v,qc{p),km)- 

Then, V = + V^'^^ + V^^^'. 
First, because 

9(^(0) +^(1)) 

d-p ^ ' 

■0'^''^ + ■0'^^' is strictly increasing in p. 

In addition, we have shown in Appendix B that 
p{i',qQ(p),k„i) is strictly decreasing in p while 
p{v,qc{p),kjn) is strictly increasing in p. Thus it is 
easy to verify that 'tp^^'>{p,X) is strictly increasing in p. 
Consequently, 4' = V'^"'' + V'^^'' is strictly increasing in 

P- _ 

b) For any A G (0, e"^], we have '0(1, A) = e^^ - A > 0, 

and 0(p, A) —5- -ho{v){k„iy~^^ < as p 0. In addition, 
the function ip is continuous and strictly monotonic in p. Thus, 
there is a unique solution p ~ i^(A) G (0, 1] for equation 

iip,~X) = 0. 

c) For given A G (0, e ^), ■0'^°) and ■0''^'* are both strictly 
increasing in p. In addition, j0'^°) = '0'-^-' = when p = 1. 
Because the solution p = a;(A) < 1, we 

V)(2)(a;(A), A) < and thus 4''^^'>{oj{X), \) = w(A)e-^(^) - 
A > 0, i.e., a;(A)e^"(^) > A. 

Appendix E 
Proof of Lemma 4 

Since the T-slot drifts of Xt can be approximated by 
the drifts of X^^^, which can be further approximated by 
closed form expressions, we start our proof by examining the 
properties of these expressions and choose the values of 7 
and S. Then by choosing T and M properly, we make the 
approximate errors small enough such that the actual drifts 
have the same properties as their approximations. 
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Notice that \ — pe^P < when p = 1 or p = /3. In addition, 
it is a continuous function and is monotonically increasing in 
p G [f3, 1]. Therefore, there exist some 5i and 7 such that 
A - pe~P < Si for all pe [/3 - 57, 1 + 67]. For A), with 
its strict monotonicity in p, we conclude that ip{p, A) < ip{l3 — 
7, A) <Oforallpe (0,^-7) and?/'(p,A) >?/'(l+7,A) >0 
for all p G (1 + 7, 00). Thus, there exists some S such that 

\- pe-P < -35/2, Vpe [/?-57, 1 + 57], (38) 
^ip,X)<-3S, Vpe(0,/3-7), (39) 

ilj{p,X)>3S, Vpe(l + 7,oo). (40) 

Now, fix 7 and 5. Using the uniform convergence of 
^d'rp{n, n, k), we can choose a sufficiently large T, such that 

for any (n, n, k) G Sjf , 



< (5. 



(41) 



where p — n/fi. 

According to Lemma 2, with the chosen T, there exists 
some M > 0, such that if n > M 01 fi > M, then 



and thus, 



^drin, n, k) - ^d'^in, n, k) 
j,dT{n,h,k) - ^d'rr{n,h,k) 



^drin, h, k) - ^JtI": ^> fc) 



<V2, 
< 5. 



(42) 
(43) 

(44) 



With ^d'rp{n, h, k) = X-pe^P, (38) and (42) together imply 



that for any (n, n, k) e 6*57, 



—dT{n,h,k) < -6, 

and thus (26) holds. 

Similarly, (27) follows by combining (39), (41), and (44); 
(28) follows by combining (40), (41), and (44). 
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